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ABSTRACT 


Shipboard measurements of the vertical profiles of mean 
humidity, temperature and wind speeds over the ocean have 
been investigated to verify their suitability to accepted 
turbulence similarity theory, In addition, the profiles have 
been used to obtain estimates of the temperature structure 
parameter. One hundred eighty-two profiles observed over 
a seven month period in 1974 aboard the R/V Acania in 
Monterey Bay have been compared to previous experimental 
results obtained in defining the meteorological parameters 
Sit ricbilonevelocity, rougnness length, drag coefficient 
and the Richardson number. The temperature structure param- 
eters computed from the profiles were provided to Lund (1975) 
for comparison with direct measurements. 

In general, favorable results were obtained in the param- 
eter comparisons which provided additional support to the 
Suitability of shipboard profile measurements defining the 
marine surface layer, However, some evidence exists in 
roughness lengths and drag coefficient values to indicate 


wooo tie ship Ss structure may be a factor influencing results. 
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I. INTRODUCTION 


The air-sea interface region has many physical properties 
yet to be fully understood, Only in the past decade has 
there been a concerted effort to make measurements in the 
planetary boundary layer over the ocean and to examine exist- 
ing laws describing it. The Navy is actively engaged in 
marine atmospheric boundary layer research for many reasons. 
One of which is to accurately access the environmental 
influence on optical wave propagation over the ocean. Exper- 
imental evidence suggests that certain near surface atmospheric 
parameters can give reasonable results in specifying the 
optical conditions for the first tens of meters above the sea. 
In this area, problems involving applied theory, shipboard 
instrumentation and data processing exist. 

A primary purpose of this work is to evaluate, for the 
marine boundary layer, expressions which have been formulated 
on the basis of previous experiments. 

Data examined in the study were measured from the Naval 
Postgraduate School research ship, the R/V Acania shown in 
figure 1. Beginning with accepted basic surface layer theory 
and advancing to turbulence similarity theory, observed data 
1S analyzed in part for its practicability in defining optical 
wave properties in an ocean environment, 

The suitability of the R/V Acania as a platform for ocean 


micrometeorological measurements was previously examined in 








Research Vessel 


Figure 1 
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detail by Cavanaugh (1974) and Welsh (1974) and verification 
of their work is a secondary objective. 

Examined in this study are ten minute time averaged 
values of wind, temperature and relative humidity measured 
at generally four levels. The observations were made during 
the months of February, March, August and September 1974 
while the R/V Acania was at anchor approximately two miles 
offshore in Monterey Bay, Coincident with these measurements 
were measurements to obtain spectra of the wind and temper- 
ature flunctuations and also measurements of optical propaga- 
tion characteristics. The agreement of the profile data with 
boundary layer theory predictions and results of others will 
be examined to determine the validity of the shipboard 
measurements, 

The boundary layer parameters examined are the roughness 
length (2), frictiomsvelecityequ,), drag coefficient (Ch), 
Richardson number (Ri), and the temperature structure param- 
eter (C7). Examinations of some of these parameters are 
accomplished by applying the data to different calculation 
methods. 

The structure parameter is calculated from an indirect 
method based on similarity turbulence theory as proposed 
by Wyngaard, Izumi, and Collins (1971). The results were 
used in comparisons by Lund (1975), who analyzed direct 


measurements of this parameter for the same data periods, 


ae 


II. THEORETICAL BACKGROUND 


A. GENERAL 

Investigations of temperature and wind fluctuations in 
the surface layer provide a basis for defining the stratum 
dynamics by turbulence similarity theory. Measurements of 
mean wind (ua), mean temperature (T) and mean relative 
humidity (q) over a time scale greater than that of the 
fluctuating motion are used to estimate the vertical fluxes 
of momentum and heat due to the fluctuations, This can be 
accomplished because profiles of the aforementioned parameters 
are maintained by the momentum and heat fluxes. Turbulence 
Similarity theory which is often referred to as the Monin- 
Obukhov theory relates the profiles to the fluctuation 
properties. Statistical analysis is used to arrive at 
empirical expressions which are accepted and are consistent. 
Dimensional analysis was used in the formulation of all 
expressions. 

The empirical relationships relating gradients of wind, 


temperature, and humidity to the two fluxes are 
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_ Z, 
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Starred quantities of wind (u,) and temperature (T,) 
include the vertical flux parameters as indicated in the 


following expressions, 


tol 


uy. = (-uw) (45), 
_ we 

ds, = - Uy (5) 

aq, = - a (6) 


B. PROFILES IDENTIFYING SURFACE LAYER CHARACTERISTICS 

Logarithmic plots of mean wind profiles can be used to 
define certain surface layer properties such as roughness 
length (20), friction velocity Cae and drag coefficient 
(C,). The latter is normally calculated for a ten meter 
height. 


Integrating equation (1) and assuming neutral or near 


neutral stability conditions, where d 2/h) = 1, yields the 
equation 
u 
— > * 2, 
Le Zz) = fs in Z (7) 


The wind speed at the surface, u(z,) is assumed to be zero 
for reasons of profile analysis. However, that may not be 
the case if the surface is also in motion, as is the case 
over the open ocean, Wind measurements taken at several 
levels and plotted as a logarithmic z profile yield an inter- 
cept of Z., and a slope of k/u,. Figure 2 is a representation 
of theoretical versus actual wind profiles in the marine 


surface boundary layer. 
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The drag coefficient is defined as the ratio of the squared 
friction velocity to the squared mean wind speed at a specific 


height, 


uy, 
oa ae (8) 
u 


C. DESCRIPTIONS OF STABILITY INFLUENCE ON THE SURFACE LAYER 
The functions Ds Pi? op are used to account for the 
stability influence in equations (1), (2), and (3). The 
stability of the surface layer regime plays an important 
role in determining flux intensities, Monin-Obukhov theory 
determines the level where mechanical and thermal production 
of kinetic energy are equal to a dimensionable length (L), 
The ratio of the height of the measurement (2) to the Monin- 
Sbuchoy length (LL) serves as a stability index, z/L. The 


Monin-Obukhov length is defined as 





2 
7 le; uy 
7 (9) 
qk Tey 
where 
ree = Ae a. Ol Ay Als (10) 


In unstable conditions the stability parameter (z/L) is 
approximately equal to the Richardson number (Ri) which is 
often used as the stability parameter, The following relation- 
Ships between the stability parameter and the Richardson 
number have been proposed by Dyer and Hicks (1970) (unstable) 


aac Webb. (1970) (stable), 
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@ is an empirically derived constant equal to ,5, 
D. INDIRECT CALCULATION OF THE TEMPERATURE STRUCTURE 

PARAMETER 

The temperature structure parameter (Cr) is the most 
important variable in determining the optical propagation 
characteristics of the atmosphere, One means for calculating 
Cr is by employing turbulence similarity theory, Since 
atmospheric refractive index fluctuations are not included 
mimenese calculations, the method indicated 1S an indirect 
estimation of atmospheric optical properties. In similarity 
theory estimates of the dissipation of turbulent kinetic 
energy (€) and dissipation of temperature variance (E9) 
provide the link between direct methods such as proposed by 
Somrsin (19591), equation 13, and the conveniently relatable 


vertical fluctuations of wind and temperature in the surface 


iavier , 


1/3 ,-5/3 


High) Baer (13) 


Research by Wyngaard, Izumi and Collins (1971) has resulted 
in nondimensionalized expressions relating atmospheric 
stability to the structure parameter. Dimensional analysis 
Pee iomNcleNie( 2) Une traction velocity (u,), the mean 


potential temperature gradient (33), Ee, €, and oe produce 
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functional expressions dependent only on the Richardson 
number. Similarity principals are elucidated by using the 
nondimensional stability parameter Richardson number defined 
as equation 14 for indirectly obtaining a comparable structure 
parameter to identify the optic properties existing in the 


surface layer. 
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Equations including the Richardson number functions 


ery , ite ay) are 
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The latter expression relates the temperature structure 
parameter and stability parameter and allows for an empirical 
Mewermination of the i3 function™at varying conditions of 
stability. If the function f3; is defined, the temverature 
structure parameter can be readily calculated from more easily 


measured parameters, a and Ri, 
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III, EXPERIMENTAL OBSERVATIONS 


eee DATA COLLECTION 

The Naval Postgraduate School R/V Acania was the obser- 
vational platform, The ship was anchored in Monterey Bay 
at a location approximated in figure 3, This location 
facilitated complimentary and simultaneous optical experi- 
ments from two shore stations, Hopkins Marine Station and 
the Holiday Inn. Normally, observational periods were three 
days in duration beginning with a one day equipment calibra- 
tion and set up, continuing with a one day on station systems 
test and finiShing with one day of data collection, Data 
collection was mainly during evening hours. The entire 
operation required considerable inter-departmental coordination 
between the meteorology and physics departments, 

Data was collected and analyzed for the dates and time 
Shown in table 1. These observations represent over ninety 
percent of all the mean profiles obtained in Monterey Bay 
during the period from 25 February to 20 September 1974, 

Figure 4 shows the general instrument configuration 
used during the experiments, Variations in the arrangement 


occurred primarily because of individual component failures, 


B. EQUIPMENT 
The mean wind measurements were made with a Thornthwaite 


Associates cup anemometer wind profile register system, 
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DATE 


26 
27 
14 
27 
28 
3 
18 
its) 
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February 
February 
March 
March 
March 
August 
september 
september 


September 


TIME 


1113-2225 
1643-2150 
1703-2138 
1823-2142 
1824-2150 
1630-2200 
1617-2202 
1724-2353 


0004-0644 


theo} ite I 
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NUMBER 
OF 


LEVELS 
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PrP Fe BP BP BP WH 


NUMBER 
OF 


PROFILES 


23 
25 
20 
ey 
14 
26 
24 
29 
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model number 104. In operation the shaft of a three cup 
anemometer unit serves as the shutter between a light source 
Momamouovoce! | for each revolution. The cups are plastic 
cones reinforced with aluminum frames, They are attached 

to the rotating shaft by stainless steel tubing spaced at 
120 degree intervals about the shaft as shown in figure 5. 
The three cup assembly sets along with the other sensors 
were positioned at four levels on the bow with electrical 
leads to the after deck house laboratory. The sets have the 
characteristics of low starting.speeds with a small amount 
of internal friction which aids in checking inertial over- 
shoot, 

Temperature sensitive quartz crystal probes, (Hewlett 
Packard model HP-2850 ) were used to meaSure mean temperatures 
at each level, RF Signals from the crystal probes and from 
a reierence oscillator were mixed in the HP-2801A readout 
unit to produce a beat freauency whose signature can be 
analyzed to within 0.001 degrees centigrade per hertz, 

Each sensor Simultaneously received pre-experiment calibra- 
tion against a platinum resistence wire thermometer in a 
Piper comeoantrolled circulating water bath over the 
expected temperature range. The accuracy in achieving a 

.005 degree centigrade correction factor was a constant for 
each probe. A 3.7 meter flexible coaxial cable is permanently 
attached to the sensor head and the mast mounted probes are 
housed in an aspirated shelter as depicted in figure 6, 


Temperature values were automatically recorded on a printer 
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tape, The Hygrodynamics Digital I hygrometer indicator and 
the Dunmore-type lithium chloride sensor were used to measure 
the relative humidity, The Dunmore-type lithium chloride 
sensor is illustrated in figure 6, The equipment operates 

on the basis of resistence change within an electrolytic 
solution generating a reference voltage variance which is 
proportional to the relative humidity change. Automatic 
temperature compensation in the instruments meet the follow- 


ing specifications for relative humidity, 


|+ 


3% relative humidity below 90% 
ap 4% relative humidity above 90% 

Calibration of the sensors was accomplished by comparative 
methods using saturated salt solutions in an enclosed 
chamber. 

Relative humidity was observed and recorded similar to 
that of temperature, as a printed output. A more detailed 
description of the devices, their calibration and deploy- 


ment can be obtained from Cavanaugh (1974), 
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IV, ANALYTICAL DATA PROCESSING 


A. DATA EDITING 

During the experiment, mean wind measurements over ten 
minute intervals were recorded in a log book, The recordings 
were the number of cup revolutions for the interval as 
recorded by electronic counters, Temperature readings were 
printed on a paper tape from a monitoring digital readout 
at approximately one second intervals, Initially relative 
humidity measurements were recorded the same as the tem- 
—Eabures. However, due to a printer malfunction, these 
readings were recorded by hand in log book after the 28 
Pepe m Cruise. Ail information was later coded and key 
punched on computer cards for further analyses, 

Preliminary screening processes were performed on 
computer listings, The data were examined for gross errors 
or inconsistencies, Very few deletions were necessary at 
this phase. Next, individual mean wind and mean temperature 
measurements were plotted to give logarithmic graphical 
profiles for each time period. The computer program devised 
for this operation also performed a least squares fit of the 
data. 

In addition, mean wind profiles for certain periods 
were closely examined for any biases which might be caused 
by the ships' hull or superstructure. Approximately ten per 


cent of all the measurements was discarded either because 
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of incongruous profiles or uncomplementary data times. 
Rejected profiles were mainly those of wind... They occurred 
principally when winds were very light. 

Temperature profiles were rejected when one or more of 
the automatically recorded readings appeared to be spurious. 


This was not discovered during the initial data inspection. 


B. DATA RECORDING ADJUSTMENTS 

Since the data was obtained over a seven month period, 
Significant alterations in data recording Soeyn place due to 
equipment failure or experiment deviations. Mean wind read- 
ings were obtained at four levels for all the periods except 
when the number of levels were two and three for March 14 
pmomeee respectively. In August, the order in which the mean 
wind levels were labeled changed from previous cruises. This 
data variation was compensated for by analyzing the measure- 
ments on a cruise to cruise basis with simple computer 
program modifications, There were no major problems encoun- 
tered in programmed analysis of the temperature data, When 
the indicated relative humidity was missing, adjustments were 
made by taking an average between the two adjacent available 


measurements, 


C. PROGRAMMING CONSIDERATIONS 

Cup anemometer readings were recorded in revolutions per 
time period, The relationship between revolutions to speeds 
in nonlinear, It changes with wind speed, A computer 


subroutine was devised to convert the number of cup rotations 
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into wind speed, The conversion accuracy for wind speeds 
decreases with decreaSing revolutions, 

An IBM System/360 source library program called LSQPL2, 
Least Squares Polynomial Fitting, was an integral part of the 
data analyses. LSQPS2 not only served to give linear func- 
tions for data profiles, but also was used as a means for 
approximating roughness length (2) values, Once the least 
Eomaces linear profiles had been calculated, profile values 
“were estimated for levels one meter either side of the 
observation levels. This was done for the purpose of defin- 
ing gradients used in computations of the Richardson numbers. 

The basis for describing stability in discerning optical 
propagation conditions is the Richardson number, In testing 
for the feasibility of applying identical relationships over 
the ocean, Richardson numbers were determined from the 
Monterey data utilizing procedures Similar to those set 
forth by Wyngaard et al, (1971), Mean profiles of temperature 
and wind were obtained by an analytical least squares fit of 
ten minute mean observations at their respective instrument 
heights. This provided the wind and potential temperature 
gradients necessary to compute a Richardson number at each of 
the four instrument levels. Equation 18 is the approximation 
for the mean wind speed gradient between two levels used in 


equation 14, 
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ou = Au = u(Z2) - u(z,) u(z (Glney) 
IZ - Az Zi ASN ZA 7S 


1. 
mpecepresents the geometric mean” height, (2, ° Ae A 
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Similar expression was used to compute temperature gradient, 
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V. DATA ANALYZATION AND COMPARISONS 


A. FRICTION VELOCITY AND ROUGHNESS LENGTH RESULTS 

Important parameters describing the marine boundary 
Mayer are the friction velocity (u,) and roughness length 
(20). The friction velocity is a measure of the momentum 
transfer described by equation 4 and the roughness length is 
a boundary parameter defined as the height where the loga- 
rithmic profile becomes a valid specification, These param- 
eters were calculated from the data by two methods. One 
method utilized all profile levels, the other utilized only 
the wind speed at one level and the air water temperature 
difference. This was done to determine what difference would 
occur if the surface layer properties were obtained using 
data from one level instead of from a four level profile. 

The first method was based on a least squares fit of a 
stability corrected four level profile. The stability 
influence was determined at each level from the Richardson 
number. The following Richardson number categories were 


used in the corrections. 


~0,05 < Ri < ,OOL neutral 
Ri > ~OoL stable 
Ri < -,05 unstable 


Stability relationships in equation il and 12 were used to 
obtain the stability length (z2/L) which was then used to 


obtain the correction to the empirical expression for mean wind, 
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The plot of the mean wind versus the logarithm of 
height with the stability correction should theoretically 
give a straight line with a slope of k/u, and an intercept 
of Zo: Equation 19 was used to incorporate the stability 


influence into the profile 


In Za-e v FD = u(z) + In Zo (19) 


k 
uy, 
The function (=) was defined as follows from Paulson (1970) 


NZ 
VF) = 2 In [Ae 5 in [Ge - 2tan ’ x +s (20) 


PG 
4 


Cl = 16(=)) unstable (21) 


Zy z 
WF) = -4.7¢ stable (22) 


A least squares fit was applied to the corrected profile to 
obtain the best estimate of the slope and intercept thus 
yielding the profile values for friction velocity and rough- 
ness length, 

milessecond method was to use agwealculation suggested by 
Cardone (1969). He formulated an iterative procedure to 
find friction velocity and roughness length given a mean 
Wind, an air-sea temperature difference and the instrument 
heights. By utilizing wind and temperature gradient expres- 
Sions combined with equations similar to 9, an iterative 
algorithm was developed, which gives stability length 
Couvergence to a desired accuracy. The friction velocity 
was obtained using the stability criterion. The roughness 


length was then calculated from the friction velocity and 
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the following expression from Cardone, 
item = C CZ) 


The constants contained in equation 23 were derived so that 
the ten meter drag coefficient (Ci) would be a minimum 
when the mean wind speed waS approximately six meters per 
second under neutral conditions, 

The results from the two methods are compared in figure 
@ which is a logarithmic plot of Z, versus u,. Comparison 


is also made with Charnock's (1955) expression, 


O a a = 81,1 (24) 

The Similar curvature of all three plots indicates agree- 
ment on the shape of the friction velocity and roughness 
ineta gwelationship, However, the least squares method 
produced roughness lengths on the order of one magnitude 
larger than Cardone's values, 

In general, both methods appear to give good results, 
however, the average values obtained by the least squares 
method show less dispersion at low friction velocities than 
does Cardone's one level method, 

A comparison of the roughness length versus the ten 
meter mean wind speed was made between results obtained by 
Ruggles (1970) and the least squares method as shown in 
figure 8. Ruggles hypothesised that the Zo peaks were 
associated with the effects of wind wave coupling. He 


designated the first peak as being coincident with the 
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minimum wind necessary for wave generation, The results of 
this study show similar maximums and minimums, most notibly 


the maximum near 3.5 meters and minimum near 4.5 meters, 


B. DRAG COEFFICIENT RESULTS 

Ten meter drag coefficients (Ci) were calculated to 
provide further verification of the compatibility of the 
experimental profile results to results of others, Friction 
velocities which were computed by the least squares method 
and also using Cardone's expression and ee ten meter mean 


wind (U5 6) were used to obtain C by the expression 


10 


u 
Cam = ae (25) 


Plots of C50 versus u, and C50 versus Richardson number 


me 
appear in figures 9 and 10, Superimposed on the graphs are 
the results of others, 

In general, some aspects of the results compare favor- 
ably with results of others and indicate that over ocean 
Broeboard app leatden of turbulence Similarity theory might 
be valid and feasible in the marine surface boundary layer. 
The least squares method produced average drag coefficients 
that were higher than the values of other experiments. This 
fomecturiabuted Co the high values of friction velocity possibly 
due to the ship's influence, The problem was addressed by 


Cavanaugh (1974), however, his results were not as pronounced 


as those in this study, 
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C, TEMPERATURE STRUCTURE PARAMETER 

The most important element for parameterizing optical 
propagation conditions in the marine surface boundary layer 
is the temperature structure parameter (Cr) The indirect 
method was adopted to define Cor at the four instrument levels 
where a computed Richardson number was available, First, it 
was necessary to use the Richardson numbers to define the 
£, function in equation 17 that had been empirically derived 
from the Kansas field data. Equation 17 was then used to 
obtain C,, 


Tomobatoreniea into the refractive idex Structure parameter, 


Subsequently, the temperature structure parameter 


Ce for additional comparisons with theoretical predictions, 
Seectral” and optical results, 
Calculated Cr 


stability conditions, Above a Richardson number of ,20, the 


decrease with height as expected for all 


f_eeunetion iam equation 17 "is undefined, 
Temperature structure parameter Ce) results were 


compared to direct measurement values by Lund (1975), 





VI. CONCLUSIONS AND RECOMMENDATIONS 


The ship's motion and Structure appeared to have had an 
effect on the wind profiles. This conclusion is based on the 
fact that the values of 25 and Cio obtained by a least squares 
fit profile data were consistently higher than values observed 
in other experiments. 

Cavanaugh (1974) observed Similar results and alluded 
toeship’sS motion influence on a logarithmic profile. Hull 
interference could disrupt the wind stream passing the 
instruments. This obstruction could cause the lower wind 
cups to slow and thus increase the apparent gradient. This 
would result in higher Zo values. This possible influence 
was examined in the analysis by disregarding the lower sensor 
level in the regression calculations. Roughness lengths were 
lowered slightly in most profiles examined, however, not 
enough to justify dropping this level in profile analyses. 

A suggestion for future experiments would be to position 
a sonic anemometer vertically over the side of the hull at 
main deck level to measure the forced vertical air motion 
caused by hull wind blockage. From data obtained in this 
manner, a Ships wind stream signature for varying wind speeds 
might be found. This could find application in the more 
complex study of measuring turbulent motion while underway. 

With regard to temperature structure parameter estimates, 


the function i is required. Sufficient data is now becoming 
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available to permit evaluation of equation 17 for over the 


2 
T 


obtained by Lund (1975) and those computed in this study 


ocean observations. Comparison of the CA spectral results 
Could aceomplish part of this task. 

Because of the quantity of data to be screened only a 
cursory analysis of the observations has been completed in 
this study. More detailed examinations into the measurements 
would certainly yield additional knowledge, particularly with 
regard to the ship's motion and structure influences and 
instrument biases. 

It would be desirable to have mean wind observations 
closer to the ocean surface. Combined with the present 
profiles, a more detailed description of wind-wave coupling 
and its effect on turbulent flow could be undertaken, The 
additional measurements could be made from a buoy positioned 
upwind of the ship. Another recommendation which is unrelated 
Pomapp lication sof turbulence Similarity theory is that a small 
portable computer be obtained for future experiments to 
facilitate reai time data analysis. The system now employed 
is extremely inefficient because many data discrepancies are 


not recognized until after the cruise, 
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